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Abstract 

Gradient theory converts the molecular theory of inhomogeneous fluid into nonlinear 
boundary value problems for density and stress distributions in fluid interfaces, contact 
line regions, nuclei and microdroplets, and other fluid microstructures. The relationship 
between the basic patterns of fluid phase behavior and the occurrence and stability of 
fluid microstructures is clearly established by the theory. All the inputs of the theory 
have molecular expressions which are computable from simple models. 

On another level, the theory becomes a phenomenological framework in which the 
equation of state of homogeneous fluid and sets of influence parameters of inhomogeneous 
fluids are the inputs and the structures, stress tensions and contact angles of menisci are 
the outputs — outputs that find applications in the science and technology of drops and 
bubbles. 


Introduction 


As witnessed by the papers presented at this colloquium, drops and bubbles, thin films 
and adsorbed layers , and contact angles are key actors in numerou ; natural and man-made 
processes. VJith our knowledge of and demands on these processes becoming more sophisti- 
cated, it is increasingly important to have a molecular level theory of structure and 
stress in interfaces. Although the formal statistical mechanical theory of inhomogeneous 
fluids at equilibrium has been developed rather extensively over the last two decades, the 
formal theory is presently intractable.^’^ Far more powerful is gradient theory,^ an 
approximation going back to Rayleigh^ and van der Waals" which was rediscovered by Cahn and 
Hilliard^ and was recently put into modern form by Eongiomo et al.^ and Yang et al.’’ 
Successes in predicting the surface tension of polymer melts,® hydrocarbons and their mix- 
tures,® and water prove that the theory is useful for real fluids. In this paper, I 

outline the elements of gradient theory and describe applications that my Minnesota 
colleagues and I have made of the theory to fluid-fluid interfaces, fluids at solid sur- 
faces, and drops and bubbles. 

Gradient theory of microstructured fluids 


A fluid microstructure is an inhomogeneous region in a fluid in which component 
densities vary appreciably over molecular distances. Any fluid is, of course, inhomo- 
geneous because of the presence of gravity. However, the inhomogeneities that result from 
gravity are so weak that component densities vary negligibly over molecular distances. 
Similarly, the inhomogeneities induced by ordinary centrifugal fields and by the temperature 
and composition gradients involved in the usual transport situations are very we.ak. If the 
component densities vary sufficiently little over molecular distances, then the thermo- 
dynamic functions can be approximated locally by the corresponding functions for homogeneous 
^luid at the local composition. In fluid microstructures the effect of the local component 
..density variations must be accounted for in the local thermodynamic functions. 

In the absence of external fields and density inhomogeneities, the Helmholtz free energy 
density is From intermolecular interactions species i and j give to f^ a contribution 

of the order of a^jn^nj , n^ and nj being component densities and a^j a characteristic energy 
parameter. The factor'^1/2 ninj is a measure of the number of interacting pairs. If the 
fluid is inhomogeneous at position r, then the number of interacting pairs in the vicinity 
of r should be corrected by some amount 1/2 5ni inj . An estimate of 6n^is rijVni, r^j the 
range of the intermolecular force between i and j and Vn^ the gradient of n^ at r. It 
follows heuristically then that the local Helmholtz free energy density of inhomogeneous 

V 1 

fluid is fQ(n) + j *'ij'^”i*'^”j terms higher order in gradients of component densi- 

ties; f_(n(r)) is the Helmholtz free energy density of homogeneous fluid at the local 
composition'n(r) t {nj^(r), n 2 (r) , • • • ,n (r) } and the terms involving density gradients are 
the Helmholtz free energy density of the' inhomogeneity. The quantity cij is prooortional to 
aijr^4, the proportionality factor arising from appropriate molecular averaging, ff an external 
conservative potential uMr) is present, then to the local free energy density is added 
t u*. Rutting together the pieces of this heuristic argument, one gets the gradient 
theoretical formula for the Helmholtz free energy F of inhomogeneous fluid: 
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c^jVn^'Vrij + ^ n^Ugld^r 


( 1 ) 


If the density gradients are macroscopic, e.g., caused by gravity or ordinary centrifugal 
fields, then Vnj^ is of the order of n^/L, L being the dimension of the system. In this case 
the quantity CijVn^'Vnj is negligible since it is of the order of (r^j/L)^ compared to the 
local value of the homogeneous fluid free energy density fn. Thus, it is appropriate to 
1 

identify .1. ^ c.^Vn^^^Vn. as the free energy density of fluid microstructure. 


Van der Waals introduced the one-component version of Equation (1) in his theory of 
liquid-vapor interfaces and Cahn and Hilliard^*** used a binary regular solution version of 
the equation in connection with interfacial structure and spinodal decomposition of sub- 
cooled homogeneous solution. In the modern statistical mechanical version of the 
theory,^' Equation (1) is derived from a formal component density expansion of the exact 
free energy of inhomogeneous fluid. Expressions are obtained which relate the local 
"influence parameters" cij of inhomogeneous fluid to the fluid radial distribution functions 
of homogeneous fluid at local component densities. The heuristic connection of the influence 
parameter c^j to a^jr^j is justified by the rigorous statistical mechanical expressions. 

In its modem version gradient theory is a very attractive description of inhomogeneous 
fluid; on the one hand the inputs fq(n) and cj^j(n) can in principle be computed from the 
molecular theory of homogeneous fluid, but on tneother hand if molecular theory is insuf- 
ficiently developed for the fluids of interest semiempirical or empirical schemes can be 
used to deduce equations of state for fo(n) and Cj^<(n). Along these lines it is encouraging 
that the molecular theoretical formulas for Cj^s ana some predictions * ^ ^ based on simple 
models imply that the influence parameters are'^often only weak functions of component 
densities. Similarly, the success of the theory with constant cj^j in predicting the sur- 
face tensions of hydrocarbon mixtures’ and water’’ argues against appreciable density 
dependence of the influence parameters. The importance of this is that one can determine 
the values of influence parameters from limited experimental data. In all the applications 
I discuss below the influence parameters are held constant. 


At equilibrium the grand potential, 

n H F - E I n^d’r , (2) 

^ V 


is a minimum in a closed system. The chemical potential, plays the 
multiplier accounting for the constraint that N^(= /nj^d’r) is fixed in a 
density distributions nj_(r), i=l,***,v, that minimize n must, according 
variations, obey the corresponding Euler equations 


role of a Lagrange 
closed system. The 
to the calculus of 


Pi = u® + p°(n) 


? V-(c..Vn^ 
J J 
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(3) 


i=l,*‘*,v. Boundary conditions appropriate for a given fluid microstructure must be 
assigned and the component densities of the microstructure determined by solving these non- 
linear differential equations. Thus, gradient theory reduces the problem of determining 
equilibrium fluid miarostruatures to a nonlinear boundary value problem. Of course, once a 
microstruccure solution has been obtained its stability has to be established by prrving 
that it is a local minimum of the grand potential SI. Typically one solves nonlinea- dif- 
ferential equations by discretization (e.g., finite difference or finite element) an^ 
iteration using the Newton-Raphson technique. ‘ ’ A biproduct of such a solution technique 
is that the matrix generated by the Newton-Raphson technique is the one required for sta- 
bility analysis of the solution so obtained, i.e., the same algorithm generates the solution 
and the elements of stability analysis. 


The pressure tensor is another quantity of interest in fluid microstructures. In homo- 
geneous fluid, the pressure is isotropic, i.e., the number of lines of force passing through 
a small element of area from molecules lying on each side of element is independent of the 
orientation of the element. This is because the molecular population is identical in all 
directions in homogeneous fluid. This is not true in inhomogeneous fluid and so the iso- 
tropic pressure of homogeneous fluid, Po(n)I, must be corrected to account for local 
component density variation. I is the unlt'tensor. Since the number of lines of force 
passing through the area element will depend on orientation, components of the pressure 
tensor are in general different in an inhomogeneous fluid, i.e., the pressure tensor | is 
anisotropic. To second order in density gradiei, .s the general formula for the pressure 
tensor is of the form^ 
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? - Po(n)J + (Vn^^) (Vn^) + + e (Vn^) • (Vtij) ]|) . (4) 

where . * * * . functions of the local composition n. The one component version of 

Equation-* (5) was^-* first proposed by Korteweg.*® The theoretical formulas for the are 

much more complicated than that for the influence parameters.^ Several simplified 
versions of the iC*-') have been investigated. ‘ ^ The simplest of these is obtained by 
assuming that the^J inf luence parameters are constant and that the n(°') must be consistent 
with the constancy. The result is^ Ij 

? “ ^o(?H ■ I iVj - J <^n^)(Vnj) + I - j(Vn^ • (Vn^)]!} . (5) 

Since fg and Pq are related by the thermodynamic relation dP^j - -d(foV) , at constant 
T, ni^*“,n^V, the gradient theory of stress given by Equation (5) requires exactly the same 
inputs as the gradient theory of the free energy. 

In what follows the influence parameters are alv;ays assumed to be constant and either 
the van der Waals (VDVJ) equation of state or one of its empirically modified successors, the 
Peng-Robinson (PR) equation,*^ is used. Both equations can be summarized as 




nkT 
1 - nb 


n‘a 


1 + i{;[2nb-(nb)-] 


( 6 ) 


i('=0 in the VDVJ equation and ij'-l in the PR equation, a and b are energy and volume 
parameters. For pure fluids the parameters are determined by the critical temperature and 
pressure for the VDW equation and these plus the acentric factor for the PR equation.*’ 

The recommended forms of a and b for mixtures is nb « t njb< and n’a •= n.n.a.i, where 


b.: and 


aj^£ are pure fluid parameters and , if* j , are^ mixture parameters to ■*bS-J determined 
by a fit of the equation to experimental dara on two-component systems. The PR equation is 
quantitatively superior to the VDVJ equation but qualitatively both are quite similar, and so 
either serves equally well the purposes of this article. 


Planar interfaces 



1 and 2 . 
ditions P, 


phases 


These component densities are of course determined by the usual equilibrium con- 
p(n(’J) = Po(o('), y°(n(*)) = u?(nV-'), i=l,*’*,v. Equation (3) can be solvec 
analytically in the case of a one-component' fluid but must be solved numerically in the 


multicomponent case.* 


The density profile of a one-component liquid-vapor interface predicted with the VDVJ 
equation at the reduced temperature kTb/a = 0.197 is shown in Figure 1. The density is 


n 

DENSITY 



Figure 1. Density in units of b” ' , distance in tc/a . kTb/a * 0.197. Ref. 14. 


^ . 
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given in units of b"‘and the distance in units of ►c/a’, a length of the order of magnitude 
of a molecular diameter. In a planar interface the normal pressure Pm, that is the com- 
ponent of pressure measured by a flat test surface lying in the intertacial plane, is 
constant in accordance with the , condition of hydrostatic equilibrium (V*P » 0, or dP{^/dx ■ 0 
since P “ P{j ii + PT.(j3 + fet) . i the unit vector along the x-axia) . On *the other hand, the 
transverse pressure^ P~, the component of pressure measured by a flat test surface orthogonal 
to the interfacial plane, is not constrained by hydrostatic equilibrium and must take on 
whatever values forced on it by the density profile. P™ cannot be constant in the inter- 
face; otherwise, the interfacial tension, which Is the aifference between the normal and 
transverse pressures integrated across the interface, namely. 




(Pjj “ P.j.)dx f 

»a> 


(7) 


would have to be zero. 

From Equation (5) it follows that 

Pt = I Po(^) + J Pn • (8) 

This result, which has been derived from several approximations’’ to the coefficients in 
the gradient theoretical pressure tensor, is heuristically very suggestive. The density 
n(x) if uniform in a plane parallel to the interface, that is the pressure is isotropic in 
such a plane. The contribution to the transverse pressure from molecules lying near the 
plane is expected to be proportional to Pg(n(x)). Molecules lying further from the plane 
are distributed to maintain a constant normal pressure Pm. The contribution of the so- 
distributed molecules to the transverse pressure should be of same order of magnitude. 

Thus, it might be argued that Pq. will be a linear combination of P^Cn) and Pjj , the 
respective coefficients being 2/3 and 1/3 reflecting the fact that the dimensions of the 
interfacial plane are 2 and the normal direction is 1 . It would be interesting to find a 
convincing derivation of Equation (8) from this point of view. An accounting of the lines 
of intermolecular force passing through a small element of area versus orientation of the 
element might be fruitful in pursuing this goal. 

An important implication of Equation (8) is that the structure of the transverse 
pressure in the interface is determined by the equation of state of homogeneous fluid. The 
normal and transverse pressures corresponding to the profile in Figure 1 are shown in 
Figure 2. As expected there is a wide region in which the interface is under tension (i.e., 
Pj < Pfj) . Since tension is positive and equals the area under the curve Pj^ - P~ versus x, 
there must be a region in which Pt ' Pm. There is however a small region of compression 
(Pt ■' Pm) on the gas side of the interface. The correspondence between the van der V!aals 
loops in the PVT phase diagram for homogeneous fluid is seen by comparing Figures 2 and 3. 
The normal pressure Pj^ is of course the liquid-vapor coexistence pressure, i.e., Pj^ = 

Po(b,i) = Pq^^b^’ ^ liquid and vapor densities, respectively. The region of 

compression in the interface arises from the region in which the pressure isotherm lies 
above the tie-lines, which locate the pressure Pj^. 

A significant feature of the theory is that the structure and stress in the interface 
are determined by the thermodynamic functions of homogeneous fluid in the metastable and 
unstable regions of the PVT diagram. It has usually been thought that in the unstable 
region of the phase diagram the thermodynamic functions are meaningless. Far from being 
meaningless, the behavior of these functions in this region is a determining factor of 
interfacial behavior. The gradient terms provide the necessary free energy to stabilize 
states chat would be unstable in homogeneous fluid.' 

As a critical point or a solution plait point is approached, the VDU loops in the 
pressure isotherm begin to flatten out and become symmetric about the tie-line, both pat- 
terns of which drive the tension cowards zero. That tension goes to zero as a critical 
point is approached is v.el 1 -known , but the mechanism of getting low tension by symmetrizing 
the VDV) loops ot the pressure Isotherm is novel. Such an example is provided in Figure U, 
in which is given density and pressure profiles of the liquid-liquid interface of carbon 
dioxide and decane. ’’ The prof iles were pre dicted with the PR equation using the mixture 
parameter values aj ^2 “ ^^11®22 ^nd Cj ^2 "0- 9>'C]^lC22 and Cj^j^ values fitted from pure fluid 

surface tension. The tension of the interface in Figure 4 is y = 0.6 dyn/cm. 
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Figure 2. Principal pressures in a planar 
interface. VDW fluid. Pressure 
in u nits of a/b‘\ distance in 
• c / a . Ref. 14. 




Figure 3. PVT phase diagram of a VDW fluid. 

Density in units of b* * tempera- 
ture in a/kb. The liquid-vapor 
densities are indicated by points 
connected with constant-pressure 
tie-lines . 



Figure 4. Density and principal pressure profiles of a C02-decane liquid-liquid 
interface of a PR fluid. kTb/a = 0.148, den.sity nj^ in units of b^‘, 
distance in «'c/a, pressure in a/b^. a, b, c CO 2 values. Ref. 19. 
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Contact anRle, wettability, and film formation 

The adsorption and wetting characteristics of fluids on solid surfaces and at fluid- 
fluid interfaces are of enormous practical interest in the design of detergents, lubricants, 
flotation and foaming agents, paints, capillary delivery devices, and the like. In spite 
of such practical importance, the molecular theory of adsorption and wetting is still a 
fledgling science, based most often on ad hoa models for various special situations. 

Gradient theory shows promise of giving a unified theoretical basis to the subject. In what 
follows, the theory is applied to adsorption and film formation at interfaces, three phase 
contact regions and the contact angle, and perfect wetting transitions. 

Suppose a fluid-fluid interface contacts a flat solid wall as indicated in Figure 5. 

If the meniscus (a mathematical surface representing the position of the interface) is 
observed at some distance R lying far enough from the solid for bulk fluid phases to exist 
on each side of the surface but not far enough for gravitational distortion to affect it, 
then the observed contact angle e obeys Young’s equation 

Y = Y „ + '( COS 0 . (9) 

' Y oc ' Y S u h 

and denote the tensions (or surface excess free energies) of fluid phases a and B 
v^ith solid phase y- Y^,c is the interfacial tension of the interface between fluids a and 
B. The basis and meaning of Equation (9), which can be derived from a force balance on the 
hemicylinder wlose cross-section is shown in dashed lines in Figure 5, has been discussed 
at length in a recent paper by Benner et al.^^ 



Figure 5. Angle of contact of fluid meniscus Figure 6. 
at a flat solid wall. 


Angles of contact of three fluid 
phases . 


If all three phase are fluid, then the menisci define three dihedral angles (Figure 6) 
obeying the force balance 


\iB _ Y B _ 

sins’ “ sins’ ' ~ sT’h0„ 

> P 


( 10 ) 


A contact angle is not always observed when three phases are brought together, 
either of the inequalities 


If 


'Yn 


' YB 


+ Y 


(B 


or 


^ ^Yt> 


( 11 ) 


then Equation (9) (or (10)) has no solution and the free energy of the system will decrease 
as the result of a thin layer of phase S (or phase a) intruding between phase y and phase a 
(or phase 6) as shown in Figure 7. The intruding phase is said to oomplctelu or perfectly 
oet the interface between the other two phases. The transition between contact angle and 
perfect wetting behavior occurs at conditions for which one of the inequalities in Equation 
(11) becomes an equality. 
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Perfect wetting Is of course essen- 
tial in applications involving spontaneous 
spreading of some fluid at an interface. 
Examples of perfect wetting are well-known. 
In the presence of air, most liquids per- 
fectly wet on clean metal surfaces , water 
on quartz, some organic liquids on water, 
some organic liquids on some polymers, 
etc. Not so well-known is critical point 
wetting, a phenomenon hypothesized recently 
by Cahn.^° He noted that according to 
critical point theory and experiment, as 
a critical point of phase a and 6 is 
approached along a temperature, pressure, 
or chemical potential path, the inter- 
facial tension approaches zero asymptoti- 
cally as 


7 


«6 



( 12 ) 


where h is the field variable <’any thermo- 
dynamic quantity being the same in all 
coexisting phases) , h^ its value at the 
critical point, and yOg a scale factor. 
Cahn postulated that thi^difference 
I'^ya ‘ will approach zero as the 

composition of the components in phases 
a and 6 approach each other, i.e., it 
obeys the scaling law 
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Figure 7. Conditions for perfect wetting 
by either phase a or phase 6. 


(13) 


sufficiently near the critical point. Since y g approaches zero faster than that postulated 
for |y - Y.,r|. Cahn concluded that there will^always exist a critical wetting value of 
the field variable, h not equal to h„, at which one of the fluids will become perfectly 
wetting. The combination of Equations (9), (12) and (13)gives asymptotic formulas for the 
contact angle at a solid as h approaches h^,^, namely. 


cose - >(_£_) 0-74 . , 

7aB 

where y^ is the value of y,^g 


^cw ~ ^c 
at h 


0.96 


cw 


(14) 


The practical significance of Cahn's theory is that one of a pair of fluids can always 
be made perfectly wetting in the presence in a third phase by adjusting field variables 
(e.g., by changing temperature or pressure or by adding some component) to get near a 
critical point. 


To test the validity of Equation (13), an unverified hypothesis, and to understand the 
relation of y and h to fluid and solid properties and interfacial structure, Teletzke 
et al.'^ studied with'" gradient theory the behavior of a one-component fluid at a flat solid 
wall. Some results of this work are of Interest here. The PR equation of state was used for 
the fluid and the wall was characterized by the wall-fluid molecule potential 


Ug(x) 


1 .d.,9 
IF 



(15) 


,i choice appropriate for walls and fluids composed of molecules interacting with the 6-12 
Lennard-Jones potential, x denotes the distance from the wall. W is a measure of the 
strength of the wall-fluid potential and d its range. 


The gradient profile equation with constant influence parameter was solved for the 
density profile n(x) by the Galerkln technique using a finite element basis set. The solid- 
fluid tension can be computed as the area derivative of the free energy, v ■ 3F/SA, or, 
equivalently, from the pressure formula, y - /'”[Pc>(nB) - Px^***’ *'b the bulk density 

or fluid far from the wall. The boundary condixions for the problem are n(x)*0 as x-*0 and 
n(x)^ng as x^»>. 

A PR fluid has an upper critical temperature Tj. equal to 0.1704 in the units a/b. 
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Gradient theory predict! a critical-wetting 
temperature Tg„ above which the liquid 
phase becomes perfectly wetting at the 
solid-vapor surface. An example is shown 
in Figure 8, for the soli d-fluid parameters 
W - 6.4 a/b and d - /c/a. Below a drop 
of liquid would not spread on the solid, 
but would form a contact angle. Above 
the drop would spread to form a perfectly 
wetting layer. The critical wetting tem- 
perature depends on W and d. At fixed d, 
decreases with Increasing W, i.e., as 
the strength of the solid-liquid potential 
increases perfect wetting occurs at lower 
temperature. At fixed W, decreases 
with increasing range d of the wall 
potential, i.e., the longer the range of 
the potential the lower the temperature 
at which perfect wetting occurs. 

Correspondin^y , the characteristic 
tension i YlV^^cw^ decreases with 
increasing ' (and, therefore, with 
decreasing U or d) . If W and/or d are 
small enough, then is near the 
critical point and Equation (14) should 
hold if Cahn's hypothesis is true. This 
was indeed found to be true (although there 
is a small detail of mean-field versus 
correct scaling laws^‘). On the other hand, 
as V and/or d increase, the perfect-wetting 
temperature T-^, decreases and Equation (14) 
no longer holds. At sufficiently low temperature the interface is very sharp (narrow) 
and the Good-Glrifalco formula*’ 



Figure 8. 


Phase and film diagram for a 
Peng-Robinson fluid. Tempera- 
ture in units a/kb, density 
in b" ‘ . Ref. 21. 


cose 



(16) 


ought to hold since it is based on a discontinuou.s interface approximation. This turns out 
to be the case. The dependence of T^^ and y^ • "^LV^^cv’^ ^ dependence of co s6 on 
Yj^y(T)/y foi various values of W are shown In Figures C and 10 for the case d - /c/a. The 



Figure 9. Variation of critical wetting 

temperature (units a/kb) critical 
wetting tension (units /ac/b*) 
with W (units a/b). Ref. 21. 



Figure 1C. Contact angle versus yv\/(.1) /y,. 
for various W. Ref. 21. 
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near-critical scaling law, Equation (14), holds for W < 3, the Good-Girifalco formula. 

Equation (16), holds for W > 9. In the intermediate range, neither formula holds. 

Several years ago, Zisman^" suggested that the contact angle correlates with the ratio 
YLv/Yc Yj, was a characteristic only of the solid, not of the fluid pK,i e. The 

implication of Figu re 9 is that this is not so unless the reduced energy par.metsr V7b/a and 
length parameter a//c7a are fixed in the series of fluids compared. As tnese ratios depend 
both on soHd and fluid, it appears that Zisman's scaling worked because the fluid parameters 
b/a and /c/a were not very ditierent for the systems compared. 

Another thing gradient theory predicts is a first order transition of the absorbed 
layer at the solid surface. In the temperature range to (Figure 8), as the bulk 
fluid density increases from zero towards the gas side of the phase diagram a composition 
is reached at which two adsorbed layers or thin-films of different thickness are predicted 
at the same equilibrium conditions. An example is shown in Figure 11. As the bulk density 
is Increased beyond the thin-film co- 
existence curve the thin-film grows 
continuously into a thick thin-film to become 
finally a perfectly wetting layer of 
liquid when the liquid-vapor phase dome is 
reached. The temperature T^.^ i^ a film 
critical point. Above as bulk density 

increases from zero to the liquid-vapor 
coexistence curve an adsorbed layer grows 
continuously through thin-film states into 
e. perfectly wetting layer of liquid phase. 

'ieiow T^u only submonolayer adsorption 
occurs with increasing density until at 
the liquid-vapor coexistence curve liquid 
appears as a drop with a contact angle. 

On the liquid side of the coexistence 
region only submonolayer adsorption occurs. 

A thin-film transition was predicted by Saan 
and Ebner^^ using an integral model free 
energy. Their thin-film coexistence curve 
is very similar to that of described here. 

The thin-film transition predicted by Cahn 
on the basis of a two-dimensional model of 
the solid is, however, qualitatively 
different. ^ ' 


It should be emphaaiaed that aaaording co the theory the patterns of film and phase 
behavior of a one-component fluid at a solid wall are general. Sufficiently near a 
critical point of a pair a and B of multicomponent coexisting fluid phases, either a or 6 
will become perfectly wetting at the interface formed by a third phase y and the nonwetting 
fluid phase. The third phase y can be solid or liquid. The critical point can be an upper 
or lower critical point. In approaching a critical point along any field variable h, e.g. 
temperature, pressure, or chemical potential of a component, there will be a critical 
wetting value h.^, a finite distance from the critical point value h^, at which one of the 
near-critical phases becomes perfectly wetting on a third phase. Outside but near tha 
coexistence composition region of a and 6. a first order tnln-film transition occurs with a 
coexistence curve lying between h-„ and h__ , h_g being the film critical point. hj.„ lies 
between h „ and hj.. If h lies between h„ and he , then sufficiently near the a-6 
coexistence region a thin-film formed between the third phase y and, say, phase a will 
thicken continuously into a layer of the perfectly wetting phase B. 

An example of the structure and stress of a liquid-vapor Interface at which a thin-film 
has almost grown into a layer of a second liquid phase is shown in Figure 12. This figure 
was taken from the work of Falls et al.^* in which theory was applied with the PR equation 
to planar interfaces and spherical drops formed in carbon dioxide and decare mixtures. The 
transverse pressure profile is highly structured, looking like that of a liquid-liquid 
interface on the left and of a liquid-vapor interface on the right. 

There is abundant evidence that the qualitative patterns of wetting transitions described 
ere are correct' ‘ ^ ^ ’ and the expected continuously thickening thin films have been 
observed by ellipsomeury . ' ' However, the predicted first order thin-film transition has not 
been verified experimentally, nor have the critical exponents of Equations (13) and (14) 
been established experimentally. 

The theory of wetting transitions and film formation at flat surfaces requires solving 
only a one-dimensional density profile problem. However, if the structure and stress of 



Figure 11. Density profile of fluid at a 
solid wall as a function of 
bulk fluid density. Ref. 21. 
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Figure 12. Density end pressure profiles of a thln-film uc a liquid-vapor 
inter face in CO 2 and decane. Density n^ ir units bl‘, distance 
in pressure in »n/b{. ComJ)onent 1 is^C02. 

kTb]_/aii - 0.148. Ref. 19. 

the three phase contact region is desired, the problem becomes of necessity greater than 
one-dimensional, even if one phase is a flat solid. If a fluid-fluid Interface contacts 
a flat solid the compon< .it densities depend on the distance x from the wall and the dis- 
tance y parallel to the wall. Benner et al," have recently studied the contact region of 
a one-component liquid-vapor interface at flat solid. The PR equation was used for the 
fluid and Equation (15) for the wall-fluid interaction potential. The space available here 
does not allow an extensive discussion of the paper. However, one Interesting feature is 
that the stress state in the liquid-vapor interface is affected relatively far from the 
Interface. An example of this is given in Figure 13. The principal pressure directions 
in the x-y plane are indicated by crosses, 
the size of which indicates the magnitudes 
of the principal pressures in the x-y plane. 

Far away from the solid and the liquid-vapor 
meniscus (defined as the position where 
n(x,y) " l/2(n_ + n^) and indicated by the 
solid curve) "the pressure components equal 
the bulk fluid value wlch is so small the 
corresponding crosses are almost invisible 
on the scale of Figure 13. At a planar 
liquid-vapor interface the normal pressure 
component (P^x) would be constant and 
therefore the pressure crosses would appear 
as horizontal lines (P 22 small) 

in the interfacial zone. Instead, even tc 
the far right of Figure 13, a distance of 
about 25 molecular diameters away from the 
wall, the normal pressure component is very 
large. Ai>d the principal stress pairs Pxx' 

P 22 undergo four sign changes in going from 
the liquid to the vapor phase. Wild and 
wonderful patterns! Are they consequential? 

A major conclusion from the work of Benner 
etal. is that at a flat surface Young's 
equation is applicable outside the contact 
region (see Figure 5) . The radius of the 
contact region is «.wo or three times the 
thickness of the interfacial zone between the 
fluid phases. of region shown is 25>c/a. PR fluid. 

Ref. 15. 



Figure 13. Principal prescare and components 
in a liquid-vapor interface near a 
wall. kTb/a - 0.1, W - 6a/b. W1 ith 
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Spherical fluid microstructures 

Visible examples of spherical microstructures are the drops and bubbles that occur in 
mists, foams, beverages, manufactured glass, and basalt lavas. Invisible examples are 
colloidal particles, vesicles, and micelles. The classical description of spherical 
structures is based on the Young-Laplace equation and the Kelvin equation expressing 
mechanical and chemical equilibrium between a bulk phase interior and a bulk phase exterior 
separated by an interface having the tension of a planar interface. In a sufficiently 
small spherical microstructure, however, the interior will not be bulk phase and the 
tension will not be that of a planar interface. Deviation from classical behavior might be 
consequential in nucleating fluids and micellar solutions, in which the equilibrium micro- 
structures are tens of angstroms in diameter and, to mention a couple of examples receiving 
special attention at this colloqulm, microdrops or microbubbles in microdrops and thin 
liquid layers on the inside of glass shells. 

Gradient theory provides a molecular level theory of spherical structures which estab- 
lishes the point at which the classical description breaks down and determines the 
structure, stress and mechanisms of stability of spherical microstructures. By way of 
example, I will outline some results of the recent investigations of Falls et ai.*'**’ on 
one and two-component mirrodrops and microbubbles. 

For spherical microstructures suspended in bulk fluid, the boundary conditions are 
3nj^/3r = 0 at r=0 and nj^(r)— >nj^n as r-*-”, n^g being the bulk phase density of i in the sus- 
pending fluid. With these boundary conditions and constant influence parameters. Falls 
et al. have solved gradient theory for a one-component VDW fluid’" and a two-component PR 
fluid.’’ Figure 14 illustrates their results for the density profiles n(r) of liquid-like 
microdrops suspended in a vapor phase. From the Young-Laplace and Kelvin equations, one 
expects the interior of the drop to be at higher pressure, and therefore higher density, 
than the saturated liquid density n®1- . This is seen to be true in Figure 14 for drops of 
radius larger than about four molecular diameters (the "radius" of the drop does not have 
a precise meaning for microdrops). However, for smaller drops, the interior density 
decreases with drop size and the density 
profile takes on a Gaussian-like shape with 
no interior bulk region. The loss of a 
bulk-like interior begins to occur when the 
radius of the drop is about equal to the 
thickness of the liquid vapor interface. 

Because the interface is curved, the 
normal pressure profile Pp, 

(P = Pjq(r)fr + P^(r) (I - ) in a 

spherical fluid structure) is not constant 
across the interface. Thus, the pressure 
profiles in a spherical interface are quite 
different from those in a planar interface 
at the same temperature (conipare Figures 15 
and 2) . This leads om to expect strong 
deviations between the tension of a 
planar interface and the tension v(R) of a 
microdrop of radius R. Similar deviations 
are expected for the Young-Laplace equation 
as classically applied, i.e. Ppj(r=0) - 
Ppj(r=»>) = 2y /R. From the thermodynamics 
or drops it follows that the appropriate 
radius R with which to describe the tension 
of the drop is the radius of the surfacs of 
tension.*’ This radius (which does not differ 
greatly from the value of r at which n(r) = 
l/2(n(r=0) + ng)) and the corresponding 
tension y(R) are predicted by gradient 
theory.’" As illustrated in Table 1, the 
tension of the drop does differ from fcr 
small drops and the classical Young-Laolace 
equation does break down. However, what is 
remarkable is that already for drops 
fifteen molecules wide the drop interface 
has virtually the same tension as a planar 
interface and the classical Young-Laplace 
equation is accurate. 
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Figure 14. Density profiles cf microspherical 
drops in a VDV’ vapor. kTb/a = 
0.197, de-’si ty in units of b" ’ , 
distance .... /c/a. 


Table 1. Microdrops in VDM vapor, 

Radius R in units of /c/a. 
kTb/a - 0.197. 



R 




2.56 

0.79 

-0.079 

2.88 

0.83 

0.017 

3.26 

0.87 

0.19 

3.76 

0.91 

0.43 

4.6-> 

0.95 

0.73 

5.52 

0.97 

0.87 

7,74 

0.99 

0.97 

10.15 

0.996 

0.99 


Figure 15. Principal pressure profiles of a 
raicrospherical drop in a VDW 
vapor. Pressure in units of a/b ~ , 
density in b“'. distance in . c/a. 
kTb/a = 0.197. Ref. 14. 

The implication of the one-component studies just outlined is that microdrop 
curvature affects the interfacial structure and the interfacial tension very little once 
the drop radius is larger than 10 molecular diameters. In multicomponent systems this con- 
clusion may or may not follow. Consider for example the bubble in liquid CO- and decane 
shown in Figure 16 (from Falls ct The radius of the bubble is only about 12 carbon 

dioxide diameters, but its component density profiles are almost identical to those of the 
planar liquid vapor interface (indicated as a bubble with R="') . 




Figure 16. Density profiles in a planar liquid-vapor interface (R=>*) 

and in a bubble (R * 12) in a CO^-decane PR fluid. Density 

n. in units of b;’, distance in »c, i /a7 , . kTbi/an = 0.148. 
R«. 19. ^ 11 11 1 11 
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The tensions of the bubble and planar interface agree within 107.. Thus, this bubble 
behaves as expected from the one-component results. 

However, as was discusseu in the previous section (Figure 12) , in a multicomponent fluid 
conditions can be such that a thin-film of an incipient third phase may be formed at an inter- 
face. These films are very sensitive to a change in a field variable. Curving an interface 
changes the chemical potential of the system (this follows from Kelvin's equation in the 
classical theory) , and so it can be anticipated that under conditions of high adsorption or 
thin-film formation the interfacial structure and stress will be very sensitive to drop 
size. Comparison of the drop component density and pressure profiles of Figure 17 (from 
Falls et u/.’’) with the corresponding planar case. Figure 12, illustrates that this is 




Figure 17. Density and pressure profiles of a drop in a C02-decane PR 
fluid. Density in units of bT', pressure in , 

distance in *c^Y^ll' = 0.148. Ref. 19. 

indeed the case. At this temperature, it has been estimated” from the theory that the 
drop structure will not begin take on the planar form until it is larger than 100 carbon 
dioxide molecules in radius. The important implication of this result is that f'n regions 
of tkin-filn formation the interfaoi a! composition oar. te greatly modified by curvature. 

This fact might be quite important in the manufacture of objects with thin, uniform layers 
of a desired material. 

Although space does not allow discussion, gradient theory predicts the work of forma- 
tion of drops and bubbles, a quantity important to the theory of homogeneous nucleation. 

In fact, the theory was used by Cahn^ ’ to support his theory of spinodal decomposition. He 
showed that the barrier expected for the nucleation of microdrops is not there owing to the 
size dependence of the interfacial tension and that as a result a homogeneous material at 
the spinodal density can transform continuously into a multiphase system. In heterogeneous 
nucleation, it is likely that thin-film formation will be an important intermediate step in 
the process when conditions are right. These matters are ripe for future work. 

Closing remarks 

In closing, 1 would like the point out a few problems to which gradient theory might 
profitably be applied next. 

1. Thin, uniform layers of fluid in a spherical solid shell. A fluid-solid potential 
will have to be introduced for shell. For uniform thin-films and thin layers of phase only 
spherically symmetric solutions need to be sought, so the problem remains one dimensional. 

By a different method Kim, Mok, and Bemat address this problem in their paper. 

2. Fluid microstructure at rough or chemically heterogeneous solid surfaces. For this 
problem a two-dimensional solid-fluid potential will have to be introduced. Thus, the prob- 
lem is two-dimensional. 

3. Contact angles with rough or chemically heterogeneous solid surfaces. This invol- 
ves the same solid-fluid potential as in Problem 2 and a fluid interface in the vicinity of 
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the solid and so is a three-dimensional problem. The problem will pose a challenge to 
computer-aided mathematics. 

4. Drop shapes on inclined rough or chemically heterogeneous surfaces. This again is 
a three-dimensional problem and an even greater challenge than Problem 3. 
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